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The explosion that results in a cosmic 7 -ray hurst (GRB) is thought to produce emission from two 
physical processes — the activity of the central engine gives rise to the high-energy emission of the 
hurst through internal shocking^ and the subsequent interaction of the flow with the external envi¬ 
ronment produces long-wavelength afterglow^^. While afterglow observations® continue to reflne our 
understanding of GRB progenitors and relativistic shocks, 7 -ray observations alone have not yielded a 
clear picture of the origin of the prompt emission® nor details of the central engine. Only one concur¬ 
rent visible-light transient has been found^ and was associated with emission from an external shock. 
Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 
minutes after the onset of GRB 041219a®. Our robotic telescope acquired 21 images during the active 
phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission as¬ 
sociated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. 
This opens a new possibility to study the central engine of GRBs with ground-based observations at 
long wavelengths. 

Prompt long-wavelength afterglow emission is predicted to arise when the reverse (external) shock 
encounters the ejecta of the explosion^’®, or through 7 -ray heating of the circumburst material^. Indeed, four 
GRBs^’^^^^ have exhibited transient optical emission that could be associated with reverse shocks, but early- 
time optical transients have not been found for the vast majority of bursts (however, new larger-aperture 
robotic optical systems have met with increasing success). Since even moderate levels of dust near the 
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GRB or along the line-of-sight in the host galaxy could effectively suppress detectable optical emission^^, 
contemporaneous observations at IR wavelengths, where light suppression is relatively minimised, offer 
a natural means to uncover any prompt emission. This was one motivation for our construction of the 
Peters Automated Infrared Imaging Telescope (PAIRITEL; figure [D- At 1.3 metres in diameter, it is one 
of the largest, completely autonomous telescope systems in the world and one of only a few capable of 
imaging at IR wavelengths (1.1-2.3 /rm; see ref. PAIRITEL acquires images with high temporal cadence 
(integration times of 7.8s) in 3 colours simultaneously. The field of view is rather large, 8.5'x 8.5', for IR 
imaging, allowing for follow-up of GRB localisations of even modest precision. 

GRB 041219a triggered the IBIS instrument on-board the INTEGRAL^^ satellite on 2004 December 
19 at 01:42:54 UTC which was reported at 01:44:05 UTC. The initial position of right ascension (RA) OOh 
24min 26s, declination (dec.) -i-62° 50' 06" was refined to 2' uncertainty at 01:47:49 UTC and a final offline 
location^ was reported at 03:31:58 UTC. The Burst Alert Telescope (BAT) on the Swift satellite^^ triggered 
and located GRB 041219a on-board at 01:42:18 UTC with a position that was within 4' of the IR source. If 
Swift had not still been in its commissioning phase with slewing disabled, the spacecraft could have slewed 
to the location within 70s of the BAT trigger. The groundbased BAT location (RA OOh 24m 37.0s, dec. -i-62d 
50'49.2") was within 48" of the IR source. As viewed by BAT, the burst duration (At) above background 
was 520s and was very bright, with up to 6.5 x 10^ cts/sec (unsaturated) between 15 and 350 keV and a 
fluence of 1.15 x 10“^ erg cm“^. The time evolution of the count rate in four BAT channels covering 15 to 
350 keV is reproduced in the lowest panel of figure|2l 

PAIRITEL began to slew on 2004 December 19 01:48:20 UTC and the first observations of the GRB 
field commenced 58 seconds later. Despite very poor observing conditions (sustained 40 mph winds, vari¬ 
able sky transmission, and 4" seeing), comparison of the first epoch of data revealed a new, variable source^^ 
not visible in the 2MASS catalogue images of the field. When compared to the astrometric grid of 2MASS 
stars in the field, we find the absolute position of the IR transient (IRT) to be RA OOh 24min 27.68s ± 0.124", 
dec. -1-62° 50' 33.501" ± 0.228", with its uncertainty dominated by the mapping to 2MASS catalogue stars. 
An optical flash was also detected in 10 images acquired before t-i-400 s by the RAPTOR experiment during 
the prompt 7 -ray emission^^. 

PAIRITEL observations of the transient continued over the following three nights, until inclement 
weather in Arizona precluded additional observations. In total, 5,790 images were acquired by the system 
over these nights. In addition, we obtained deep J-band imaging on December 20 and 21 UTC using the 
NIRC-I^^ instrument on the Keck I 10 metre telescope on Mauna Kea, Hawaii. Consistent within the 
astrometric accuracy of the IRT from December 19, we found a point-like source in J-band (figureQ. Over 
these two nights, that source was seen to fade by 1.0 magnitudes, confirming its identification with the IRT. 
The Keck images revealed two sources within 2.5" of the transient position (SI: J ^ 19.7 mag, 2.5" North- 
north-east; S2: J 21.4mag, 1.5" East): both were unresolved apparent point sources. The source SI is 
bright enough to contaminate the PAIRITEL J-band aperture photometry on Dec 21, which accounts for the 
difference between our measurements and the fainter measurements from Apache Point^° and Keck on the 
same date. When comparing PAIRITEL photometry with higher-resolution Keck and APO results for Dec 
21, the flux from S1-I-S2 appears to be a 51% contamination in J7-band, 58% in J band, and has negligible 
contribution in the Kg band. 

The resulting light curves (see Supplemental Table #1) shown in figure |3] reveal a complex time 
history of the afterglow. The first six PAIRITEL exposures at f-i-7.2 min after the trigger show a source that 
brightens, then fades very rapidly in all filters by about t+9 min, and then rebrightens by f-i-20 min. Using 
our data and data reported in the literature, we fitted the light curves as the sum of three smoothly-connected 


2 


rise and fall brightening events; the results of these fits are shown in figure |3] During the first few days, 
the source colours, though rather uncertain, appear consistent with a single value of the spectral slope of 
(3 ^ 0.4 (figure 4). Additionally, there is some evidence that the IRT was redder during the “flash” event at 
t+1.2 min. 

How might the light curve be understood as emission from the reverse and forward shock? The 
electrons in the shock are assumed to be accelerated to a power-law spectrum with number density as a 
function of energy {E) proportional to E~'p. In a constant-density circumburst environment, a reverse shock 
is expected (see ref. to rise rapidly (a = 3p — 3/2, with f^, oc t°‘) and then, in the “thin shell” case 
(see below), decline with a = —{27p -|- 7)/35 Ri —2 after the emission peak, corresponding to the time 
(tx) that the reverse shock crosses the explosion ejecta. The IR emission from the forward shock^^ should 
rise slowly as a = 0.5, then decline as a = 3(1 — p)/4. Associating peak 2 with the reverse shock and 
peak 3 with the forward shock, we hnd reasonable agreement, to within the measurement uncertainties, of 
the data with this model. In particular, the three implied values for p (2.5 ± 1.1, 4.2 ± 4.1, 2.6 ± 0.2, 
for the reverse-rise, reverse-decline, and forward-decline, respectively) are all consistent with the usual 
range^ of p = 2.2 — 2.5. A consequence of this interpretation is that, in the absence of effects due to 
collimation of the burst, radio afterglow emission should be dominated by a rising reverse shock, peaking at 
time f = (10.2±1.8)(i//8.4 day; this is thus far conhrmed with reports of rising radio emission 

at least to day 2.9 (see ^^). 

With this interpretation, there are three puzzles. First, we would expect^ the source to have become 
bluer during the forward shock rise, which is not required (though is not excluded) by our data. The second 
puzzle concerns the timing of the reverse shock peak relative to the GRB duration. We expect^^’^^ fx = 
1670s > At only when At is less than the time when the shock begins to decelerate (commonly deemed 
the “thin shell” case). This deceleration time occurs when the shock has swept up from the circumburst 
environment a quantity I/Tq times the entrained mass, where Tq is the terminal Lorentz factor of the shock. 
A delayed reverse shock crossing time requires Tq < 67£'52*n“^/®(Af/520s)“^/® [0.5(1 -|- z)]^/®, where n 
is the circumburst particle density in units of baryons cm“^, E ^2 is the energy in the shock in units of 10 ^^ 
erg, and the redshift of the burst is z. This initial Lorentz factor constraint is uncomfortably smaller than 
the limits placed on previous values^^ of Tq, which suggests that either the burst occurred at high redshift 
(which is excluded by the RAPTOR detection^*), was exceptionally energetic, or occurred in a low-density 
environment. Alternatively, the observed variability could be due to inhomogeneous density structure in the 
circumburst environment or delayed energy injection into the blastwave. 

The third and most intriguing puzzle is the physical origin of the IR flash at fl = 462s that has marginal 
evidence for appearing more red than the afterglow at later times. Regardless of the interpretation of events 
2 and 3, the rapid rise and fall appear to preclude an association with a reverse shock: the fit of a power-law 
decay slope to the Kg data acquired less than 10 min from burst trigger yields a = —18 ± 5 whereas setting 
a = — 2, as expected of a reverse shock, yields an unacceptable fit (reduced = 4.4). The duration of the 
IR flash {6t « 45s, taken as the full-width at half-maximum of the model fit to the data) is comparable to 
the widths of the largest timescale for substructure in 7 -rays; indeed, this 6t is remarkably similar to that of 
the optical pulse^ in GRB 990123. The ratio of the width to the time after trigger t ~ 462s, St/t = 0.10, 
is similar to that seen in individual pulses in bright GRBs^^ but a factor of ~10 smaller than the St/t of the 
optical flash of GRB 990123. 

We suggest, therefore, that the origin of the first peak is from the internal shock that itself produced 
the GRB. Such emission is possible if the synchrotron cooling frequency from the internal shock emission 
is well below 7 -ray frequencies^. Indeed, if the IR emission is due to internal shocks, then the observed 
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flash may be due to a superposition of several unresolved shorter-timescale pulses. Future observations of 
IR flashes in the Swift era will no doubt test the ubiquity and nature of rapidly variable early-time emission 
as presented herein. 
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Figure 1: Images of the IR flash associated with GRB 041219a. The discovery was made with PAIRITEL, located 
atop Mt. Hopkins, Arizona. Normally, the telescope’s observing plan is autonomously scheduled before nightfall by a 
routine that optimises the preset priorities and scheduling constraints of all objects in the PAIRITEL database. When 
the new GRB alert was received, a series of observations were automatically inserted into the observing queue and the 
telescope began to slew to the target held. Left: Ealse-colour image of the IR Elash (yellow arrow) of GRB 041219a 
inside the 2' (90% confidence) INTEGRAL error region* (yellow circle). At early times, the source is the reddest 
object in the held, indicative of very high extinction due to dust in the disk of our Galaxy. Middle six panels: The 
Two Micron All Sky Survey (2MASS) images of the held from 2000 June 15 UTC compared with the three colour 
images showing the infrared transient (IRT) several minutes after the GRB triggered. Time relative to the GRB trigger 
is given as is the number (N) of individual images combined. The individual 2MASS images are 1.3 s integrations, 
so the combination of 6 2MASS images is equal to a single PAIRITEL image. The 2MASS and PAIRITEL fields are 
centred on the IRT and are approximately 1' on a side (shown as a blue box in the colour image). A local background 
from a 2-D median was subtracted from all images to remove large-scale background variations. Eor all images North 
is up and East is to the left. Right: The Keck NIRC-1 imaging at J-band on December 21 UTC, in the same region 
as the PAIRITEL images. NIRC-1 images were combined in the usual manner and the combined images had a seeing 
of less than 0.8" on both nights. The IRT, as well as two unresolved nearby sources (SI, S2), are labelled. 

High resolution version of this figure: http : / /astro . berkeley . edu/~ jbloom/grb041219a/figl.ps 
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Figure 2: Evolution of the IR flash associated with GRB 041219a. We argue that the initial pulse is not traditional 
afterglow from an external shock but instead related to the internal shocks of the central engine. Shown are measure¬ 
ments derived from individual 7.8 s exposures as small, light-coloured points, and detections from stacks of images as 
heavy, dark-coloured points; error bars are 1 a, estimated from photon noise and the distribution of randomly-placed 
apertures on the images. Also plotted is the light curve from the 15-350 keV channels of the BAT instrument aboard 
Swift. Reductions: The response to variable sky and bias in the detectors was estimated for each exposure in each 
band by median-combining all the exposures taken within a two to four minute window. A flat-field correction for fixed 
pixel-to-pixel variations in detector gain was made using images of the dawn sky. Subtracting the dark-rsky response 
and normalising by the flat-field produced a reduced image. Photometry: Reduced images were measured either 
individually or as stacks of co-added images. Since individual images undersample the seeing due to atmospheric 
blurring, photometry was performed differentially in an aperture of fixed size. All of the images were aligned to a 
common reference image to an accuracy of approximately 0.1 pixel. Images were measured individually or in stacks 
created by summing individual images with weights determined by the signal-to-noise of each image. The measured 
flux at the position of the GRB was compared to the flux measured for a set of nearby comparison stars. Magnitudes of 
the comparison stars are known to a high accuracy (2%) from the 2MASS catalog. PAIRITEL uses the same detectors 
and filters as 2MASS, so our differential photometry is expected to be free from systematic offsets. 
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Figure 3: Long-timescale variability of the IR transient associated with GRB 041219a. The iTg-band light-curve is 
reasonably fit (grey curve) as a sum of three power-law rise and fall events with peaks at times since the burst [fl, t2, tS] 
= [7.73m ± 0.14m, 27.5m ±4.9m, w500m] and specific brightness fi,(t2), /j/(f3)] = [1.6 ± 0.3 mJy, 2.5 ± 

2.8 mJy, 0.5 ±0.1 mJy]. The rise and fall near the second peak is fitted by a 2 ,r = 6.1 ± 2.9 (a 2 j = —3.4 ± 2.8), 
with the specific brightness changing as f^, ex The final rise and fall parameters depend on the poorly constrained 
value of fa. Fixing = 500m, we find aa = 0.3 ± 0.1 and asj = —1.2 ± 0.1. The first peak has both rapid 
rise and fall times (see text). The global per degree of freedom (dof) is an acceptable 1.24. The reported 1-a 
parameter uncertainties do not reflect the covariance between the parameter fits. Additionally, all of these derived 
values depend on the choice of a smoothing parameter to connect the two power-laws (s in eq. 2 of Nakar & Piran^^; 
here, we chose s = 5). The J- and iJ-band light curves are consistent with the ATs-band shape, aside from the overall 
flux normalisation. Photometry; Large numbers of individual exposures were combined with sub-pixel sampling and 
produced a high signal-to-noise, high-resolution image for each epoch. The point-spread function in these images is 
well-sampled, so the technique of difference image photometry was utilised. The uncertainty in the differential fluxes 
of the IRT between individual epochs is near photon-limited. The differential fluxes are converted to magnitudes from 
aperture photometry at the position of the IRT in the late-time (Dec 21) reference image. Other data, plotted as open 
points, include two J-band observations made with the Keck I telescope, JHKg observations from Apache Point 
Observatory^®, and three Kg measurements (see and re^rences therein) from Mt. Palomar. 









Figure 4: Colour evolution of the afterglow of GRB 041219a, corrected for extinction due to interstellar dust. The 
dust extinction^^ at the location of GRB 041219a is high A(V band) = 5.9 (and highly uncertain). Using a common 
dust extinction law (Ry = 3.1), the extinction in magnitudes is A{J) = 1.6, A{H) = 1.04, and A{Ks) = 0.66. We 
fitted (solid squares) for a spectral index /? with our J, H, and Kg measurements at each epoch using the flux relation 
F„{t) = Fo{t){\cl\ hY, where Ac is the central wavelength of the filter bandpasses (c = [J, FI, K]). Including 
the 2 :-band detection^^ (Ac ~ 9100 A) that was nearly coeval with our IR observations, results in a steeper spectral 
slope of 0.80 ± 0.39 (open square) than found with the IR data alone. Our fit to the Apache Point JHK data^*^ is 
also shown with an open triangle. With fits that yielded y^ldof greater than unity, we scale the resulting 1 a errors 
by Though there is some evidence for secular trends, the resulting values for (3 are consistent with a single 

value of /3 w 0.4. We caution, however, that this common value of the spectral index is strongly dependent upon the 
assumed dust column. For different values of the extinction, A{V) = 4.0 and A{V) = 7.0, the fits to the spectral index 
are (3 w 0.7 and (3 w 0.2, respectively. Indeed, the standard forward-shock model^^ tends to favour a larger value 
(/? w 0.55) than found with the IR data alone. For the first peak (during the GRB), there is marginal evidence that the 
source is intrinsically redder than at later times. 
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